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ABSTRA  CT.  The  aim  of  this  paper  is  to  present  an  approach, 
able  to  deal  with  all  possible  connections  of  voltage  and 
current  sources  and  impedances,  combining  conductors  in 
which  the  skin  effect  is  taken  into  account  and  conductors  in 
which  skin  effect  is  neglected This  approach  is  obtained 
using  a bottom-up  methodology.  In  this  way,  the  meaning  of 
terms  in  the  generalized  approach  is  naturally  inherited  from 
some  specific  problems.  This  model  is  presented  in  a compact 
form,  preserving  sparse,  symmetric  and  positive-definiteness 
matrices . The  vectors  and  matrices,  computed  during  the 
solution  stage,  are  employed  in  engineering  calculations  as 
current  and  inductance  computations  providing  compact 
expressions  suitable  for  efficient  algorithms.  Finally,  the 
proposed  approach,  implemented  in  a F EM  software  package 
developed  by  the  authors,  is  applied  to  the  study  of  a three- 
phase  transformer,  supplied  with  a balanced  three-phase 
voltage  (sinusoidal  and  nonsinusoidal)  and  loaded  with  an 
unbalanced  three-phase  RL  impedance.  The  agreement 
between  the  computed  and  experimental  results  shows  the 
validity  of  the  proposed  model  and  its  implementation. 

1 INTRODUCTION 

The  natural  excitation  for  the  Finite  Element  Method  (FEM) 
in  electromagnetism  is  current  (current  density).  This  current 
is  supplied  by  an  external  electric  circuit;  and,  usually,  the 
value  of  this  current  is  unknown,  as  the  circuits  consist  of 
passive  elements  and  voltage  sources.  Therefore,  FEM  and 
electrical  equations  are  coupled  by  these  unknown  currents. 
There  are  two  different  approaches  to  couple  these  equations. 
In  the  first  approach,  called  direct  coupling,  the  equations 
obtained  from  FEM  and  the  electrical  circuits  are  integrated  in 
a single  system  of  equations  having  the  magnetic  vector 
potential  as  its  solution.  In  the  second  one,  called  indirect 
coupling,  two  different  systems  of  equations  are  obtained;  one 
of  them  takes  into  account  the  FEM  equations,  and  the  other 
one,  the  circuit  equations.  In  this  paper,  a method,  based  on 
the  direct  coupling  of  FEM  and  circuit  equations  is  presented. 

Many  approaches  on  coupled  2D  Finite  Element  models  and 
circuit  equations  have  been  proposed  f l]-[7].  Most  of  them, 
take  into  account  the  skin  effect  for  all  the  conductors  or,  they 
neglect  its  influence,  without  evaluating  its  magnitude  over 
each  conductor.  In  this  paper,  a bottom-up  methodology  is 


applied,  generalizing  the  meaning  of  the  studied  examples, 
instead  of  a theoretically  obtained  approach.  So,  some 
common  winding  connections  were  studied— one-phase  circuit 
and  three-phase  circuits  with  star  and  delta  connections—, 
being  presented  here  the  three-phase  star  connection  handling 
unbalanced  load  and  voltage  conditions.  In  this  paper,  the 
model  is  applied  to  the  study  of  a power  transformer,  although 
it  has  been  applied  to  other  non-static  machines  as  induction 
motors,  [10]. 


2 BASIC  EQUATIONS 

Assuming  that  the  displacement  currents  are  neglected, 
Maxwell’s  equations  using  a magnetic  vector  potential 
formulation  can  be  written  in  2D  as: 


v(B) 


SA(x.yA 

8y  J 


-J(x,y)  (1) 


Where  A(x,y)  is  the  magnetic  vector  potential,  J(x,y)  the 
current  density  (z  component)  and  v the  reluctivity  that 
depends  on  the  magnetic  flux  density  vector  (B).  In 
anisotropic  materials,  reluctivity  depends  on  |£|  and  the  angle 
(0)  between  B and  the  rolling  direction  of  the  material 
v=v(^4,0). 


The  expression  of  current  density  depends  on  how  the 
conductors  are  defined.  So,  the  conductors  can  be  modelled  as 
solid  conductors,  in  which  case  the  skin  effect  is  taken  into 
account,  or  stranded  conductors,  in  which  case  the  skin  effect 
is  neglected.  Stranded  conductors  are  particularly  useful  to 
model  windings  in  which  an  individual  mesh  considering  their 
relations  would  exceed  in  many  cases  the  requirements  of 
actual  computers  without  providing  a significantly  greater 
accuracy.  If  all  the  conductors  in  the  problem  are  solid 
conductors,  the  following  equation  is  employed,  [8]  [9]: 

r)  A 

V X vV  X A = crEeit  - cr  — - ovc  x V x A (2) 

St 

The  first  term  on  the  right  side  in  (2),  Eext  being  the  external 
electric  field  and  g the  electrical  conductivity,  represents  the 
current  density  due  to  the  voltage  source.  The  second  term 
takes  into  account  the  current  density  due  to  the  temporal 
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variation  of  A.  The  third  one  represents  the  current  density 
due  to  the  movement,  being  ve  the  velocity  of  the  conductor 
with  respect  to  B.  By  using  a frame  of  reference,  such  that  the 
relative  velocity  becomes  zero,  (2)  can  be  written  as,  [8]: 

d\ 

V x vV  x A = <rE  t - <j — (3) 

dt 


The  total  current  of  a conductor  can  be  obtained  by  integrating 
the  right  side  of  (3): 


In  2D  problems  (3)  becomes: 


Figure  1.  The  windings  are  star  connected. 


-V(vVA)  = cr--  <j  — 
l dt 


(5) 


Where  ux  is  the  voltage  between  the  terminals  of  the  conductor 
and  / the  length  of  the  conductor  (thickness  of  the  considered 
domain  along  the  z-axis).  Equation  (5)  takes  into  account  the 
skin  effect. 

When  the  conductors  are  considered  stranded,  it  is  supposed 
that  the  current  (if)  is  the  same  in  all  the  conductors  connected 
in  series  (winding);  then,  their  current  density  can  be  obtained 
by  J=Nif/S.  Where  S is  the  entire  section  of  all  the  conductors 
(N turns)  integrating  the  winding,  o the  electrical  conductivity7 
and  Rf  the  DC  resistance.  The  current  can  be  calculated  in  2D 
problems,  by  considering  the  voltage  K/r  between  the  terminals 
of  the  winding  as  the  sum  of  the  voltages  of  its  conductors: 


7 


lN  jf^-dS 


Rf  R,  S J ■ ’ dt 


Using  (1)  and  (6): 

Nuf  aN_ 

SRf  ' S 


-V(  yVA)  = 


IT  — dS 

JJ  dt 


(6) 


(7) 


Summarizing  (3)  and  (7)  into  a single  equation  for  2D 
problems: 


NUr 


SR 


f 


crN 

S 


s 


dA 

dt 


us  dA 

cr  — -<r  — 
/ dt 


dS  Stranded 


Solid 


(8) 


3 EXAMPLE.  THE  WINDINGS  ARE  STAR 

CONNECTED 


connected  (figure  1).  They  are  modelled  with  FEM,  and 
indeed  can  or  can  not  share  the  same  2D  magnetic  circuit.  The 
total  DC  resistances,  including  the  resistances  corresponding 
to  the  windings  and  the  load,  are  called  Rr,  Rs,  Rt.  In  the 
same  way  LR,  Ls,  LT  take  into  account  the  inductance 
corresponding  to  the  load  and  the  part  of  leakage  flux  not 
considered  in  the  two-dimensional  model. 


From  figure  1 the  following  equations  can  be  obtained: 


Rrir  + Lr  - [ys  + Ls 

^ )'(^/r+Ar  *■) 


VST  es  eT  + 1 Rsis  + Ls 


iR  + is  + ir  - 0 

values  of  eR , es  and  eT  can  be  obtained  using: 

rdA^ 


er-i  NJ 


dt 


j = R,S,T 


(9) 

(10) 

(11) 

(12) 

(13) 


where  Am  is  the  mean  magnetic  vector  potential  on j winding, 
defined  by: 


fffij  6K2 


(14) 


•{'*(')}-(»)' MW) 


The  authors  have  solved  some  complex  winding  connections. 
In  this  case,  three  windings  of  stranded  conductors  are  star 
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Being  {w}T  the  phase  vector,  which  provides  information 
about  the  topology  of  each  winding,  (p(x,y)  the  shape  or 
interpolation  functions,  whose  formulas  are  well-known  for 
classical  elements,  and  Q the  entire  domain  where  the  integral 
is  extended  by  using  the  phase  functions  fiy. 


l 

-l 

o 


V (x9y)  € 5; 
If  C x9y ) e Sj 
Otherwise 


(15) 


Using  the  Crank-Nicholson  time  stepping  method,  the 
following  equations  are  obtained: 


it+At 


At 


(16) 


1 fdA 


21  dt 


+ 1- 


dK 

dt 


■A' 


At 


(17) 


Using  (9),  (10),  (1 1),  (12),  (13),  (16)  and  (17): 


.t+At  1 

jr  ; 

_ 

Mr  0 ~Ms 

\ V'rs^  + Vw 

.t+At  \ 

IS  1 

~Mt  Mr  0 

| rSf+VsT 

.t+At 

It  J 

Mi 

0 - Mr  Ms_ 

[vif+Vn 

(18) 


Ms  + Mt 

• Mt 

- Ms 

~Nr 

0 

Mt 

Mt  "t~  Mr 

Mr 

0 

Ns 

Ms 

Mr 

Ms  + M r_ 

0 

0 

jt+At  At  ) 

Am.R  ~ Am.R 
AZs  ~ A*mJ5  [ 

AZt-ALt  J 


‘ M'r(Ms  + Mt ) 

M ' sMt 

“M  'tMs 

M'rMt 

M 's(  Mt  + Mr) 

M 'tMr 

M'rMs 

M'sMr 

M't(  Ms  + Mr) \ 

l 

Using  (20),  (21)  and  the  definition  ofAm  in  (14)  the  following 
equation  is  obtained: 

-WIWM 

Where  the  term  [ W]  is  called  the  phase  matrix  defined  as  the 
matrix  containing  the  phase  vectors  presented  in  the  problem 
(three  in  this  case). 

Using  (1)  and  (22): 

([s'^]+/[^]I^][Gi][Ar][»']r){r4'}= 

(23) 

where  [St+At]  is  the  stiffness  matrix.  Therefore,  equation  (23) 
models  2D  magnetic  circuits  with  a three  star  connected 
windings  of  stranded  conductors,  the  magnetic  vector 
potential  in  the  nodes  at  t+At  being  the  unknown  values. 
Besides,  current  in  the  windings  at  / must  be  computed  using 
efficient  algorithms. 

4 THE  GENERALIZED  SYSTEM  OF 
EQUATIONS 

The  example  analyzed  above  can  be  generalized  to  different 
types  of  connections  and  conductors  (solid  and  stranded)  [1 1]. 


where: 

M,  = RiAt  + 2Li  j = R,S,T 

,1  j j j 7 7 

M',  = RlAt-2Li  j = R,S,T  (19) 

J J J J 7 

Ml  = MrMs  + MsMt  + MtMr 
Equation  (18)  can  be  written  in  a packed  form  as: 


^]+^[r]+/[^][Al][G43[Al][irfj{^}= 

(24) 

+L]hc]K:T}+^[^][^][Gj{«rA')+'} 

-([^][^][gJ[gJ-[iu][^]){;l} 


A/[Gj{«rw} 

The  current  density  can  be  obtained  using: 


X 

0 

0 

Sr 

Ns 

c 

,/^A/  J 
IR 

fs* 

> = 

0 

0 

is 

Tt+tJ 
. •*T 

0 

*4 s 

0 

1 

.t+At 
{ lT 

So,  by  example,  if  a specific  problem  includes  windings  of 

(20)  stranded  conductors  with  star  connection  and  other  windings 
with  whatever  connections,  the  contribution  of  the  star 
connected  windings  to  the  terms  of  (24)  is  obtained  from  (23). 

The  matrices  in  (24)  maintain  the  sparsity,  symmetry  and 
positive-definiteness  of  the  traditional  FEM.  Obviously, 
sparsity  depends  on  geometry  and  the  number  of  nodes  since 
the  mean  magnetic  vector  potential  introduces  connections 

(21)  among  the  nodes  of  the  windings  modelled  with  stranded 
conductors  (14). 

Defining  «,  m and  p as  the  number  of  nodes,  and  windings 
made  of  stranded  and  solid  conductors  respectively,  the 
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meaning  and  dimensions  of  the  terms  in  (24)  can  be  expressed 
as: 


• [‘SJnxn  is  the  stiffness  matrix. 

• [^Inxn  is  called  the  mass  matrix,  containing  information 
about  the  current  distribution. 

• [^Inxm  is  the  phase  matrix  showing  information  about  the 
distribution  of  windings  made  of  stranded  conductors.  It  is 
composed  of  the  phase  vectors  defined  in  (14). 

• [Ga]mxm,  [Gb]mxm  and  [Gc]mXm  take  into  account  how  the 
windings  defined  in  [W]  are  connected.  In  (18)  and  (20) 
they  are  defined  to  star  connections. 

• {At+Ai}  n is  the  vector  of  unknown  nodal  values  representing 
the  value  of  the  magnetic  vector  potential  at  the  current 
time  step  t+At 

• {A'} n represents  the  value  of  the  magnetic  vector  potential 
at  the  previous  time  step  t. 

• [Q]  has  the  same  meaning  in  solid  conductors  as  [W]  in 
stranded  conductors. 

• ] represents  the  addition  of  voltages  applied  to 

the  solid  conductors  at  the  current  time  step  t+At  and  at  the 
previous  time  step  t. 

• Jm  represents  the  addition  of  voltages  applied  to 

the  stranded  conductors  at  the  current  time  step  t+At  and  at 
the  previous  time  step  t . 

• {icab}m  represents  the  value  of  current  at  the  previous  time 
step  in  conductors  defined  as  stranded. 

5 COMPUTATION  OF  ENGINEERING 
RESULTS 

As  shown  in  (23)  and  (24),  currents  in  windings  must  be 
calculated  during  the  solution  stage.  Besides,  this 
computational  cost  must  be  kept  as  low  as  possible.  Other 
engineering  results  such  as  inductance,  magnetic  flux  and 
electromotive  force  can  be  obtained  by  making  use  of  the 
terms  defined  in  (24),  computed  at  each  time  step  avoiding 
time  consuming  calculations  in  the  postprocessing  stage. 


Current  in  windings  made  of  conductors  considered  as 
stranded  ones  can  be  calculated  using  (1),  (5)  and  the 
definition  of  phase  vector: 


[5'+A']{^'+A'}  = {/'+A'}  = X^{<yt}  ;/+4'  (27) 

bl 


Where  m is  the  number  of  windings.  The  current  of  a winding 
(/)  can  be  computed  multiplying  both  terms  in  (27)  by  the 
transposed  i-phase  vector  {(Oj}T,  taking  into  account  that  the 
product  of  {o)i}T  and  {ga}  is  null  when  k is  not  the  same  as  i: 


N‘  mr  m 


(28) 


This  method  computes  the  phase  current  using  as  data,  the 
phase  vector,  the  number  of  turns,  together  with  the  stiffness 
matrix  and  the  vector  of  the  magnetic  vector  potential  nodal 
values  at  the  current  time  step/+^k.  Therefore,  all  the  terms  in 
(28)  are  computed  during  the  solution  stage  of  FEM. 


5.2  Flux  and  Electromotive  Force  Calculations 

These  results  can  be  easily  computed  using  (13)  and  (14). 


5.3  Inductance  Calculations 


Inductance  has  a geometrical  meaning  when  all  the  materials 
exhibit  a linear  behaviour.  In  other  case,  several  non- 
coincident definitions  can  be  taken  into  account  [12].  Using 
the  definition  based  on  flux  linkages,  self-inductance  can  be 
calculated  in  windings  made  of  stranded  conductors  as: 


lN,{oX\ 

f J'+A'l 

r j 

1 

H1  [>*; 

(. A '+*} 

NM)1  {<»,} 


=/^{®T[s'+aT>} 


5.1  Current  Calculations 

The  winding  currents  can  be  calculated  using: 

if(t)=  Jj(t)  dQ 

tv 


In  solid  conductors  (25)  can  be  expressed  as: 


In  the  same  way  mutual  inductance  can  be  computed  using  the 
flux  linkage  definition: 


(25) 


Mu 


MXX-  inMAa"*} 


. t+& 


N,{wi\  { co, } 

= / 1 W 


(30) 


(26)  As  the  current  is  computed  in  (28),  self  and  mutual 
inductances  are  calculated  using  the  phase  vectors  and  the 
stiffness  matrix,  that  can  be  dependent  on  time.  Stiffness 
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matrix  is  affected  by  the  B-H  curve  when  it  exhibits  a 
nonlinear  behaviour. 

In  [13]  the  inductance  computation  in  case  of  solid  conductors 
is  stated. 

6 EXAMPLES 

A software  package  has  been  developed  implementing  (24).  It 
has  been  tested  with  many  examples  of  increasing  complexity 
over  transformers,  considering  several  balanced  and 
unbalanced  connections  in  primary  and  secondary  -one  phase 
and  three  phase  transformers-,  steady  and  unsteady  states, 
sinusoidal  and  nonsinusoidal  voltages  and  current  sources 
(balanced  and  unbalanced);  two  of  these  examples  are 
presented  below. 

6.1  Steady-state 

The  model  is  applied  to  a 3 kVA  three-phase  AY  transformer 
(figure  2).  The  transformer  is  supplied  with  a balanced 


Load 


Figure  2.  The  transformer  is  delta-star  connected. 


0 0.005  0.01  0.015  0.02  0.025  0 03 

Time  (s) 


Figure  3.  Line-to-line  supplied  voltages. 


nonsinusoidal  three-phase  voltage  (figure  3),  generated  with  a 
programmable  AC  power  source.  These  waveforms  are 
applied  to  the  software  package  using  files  obtained  from  an 
oscilloscope.  The  secondary  is  loaded  with  an  unbalanced 
three-phase  RL  impedance.  Figure  4 shows  the  simulated 
currents  in  steady  state,  and  figure  5 shows  the  computed  and 
measured  currents. 


Figure  4.  Primary  and  secondary  computed  currents. 


14  15  16 


Figure  5.  Computed  (solid)  and  measured  (dashed)  currents. 
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LAB/Primary  line  currents  - Time  (s) 

Figure  6.  Computed  (FEM)  and  measured  (LAB)  currents. 


6.2  Non  steady-state 

The  transformer  is  now  YY  connected,  and  loaded  with  the 
same  unbalanced  three-phase  RL  impedance.  Figure  6 shows 
the  computed  and  measured  primary  currents;  the 
discrepancies  between  them  are  due  to  the  unfavourable 
conditions  of  the  high  supply  voltage  and  the  instant  chosen 
for  the  switching-on  of  the  transformer. 

7 CONCLUSIONS 

A generalized  approach  for  FEM  and  external  circuits 
coupling  has  been  presented  considering  both  types  of 
conductors  simultaneously,  stranded  and  solid  respectively. 
This  model  has  been  established  using  a bottom-up 
methodology.  The  meanings  of  the  terms  of  the  particularized 
equation  (23)  have  been  generalized  to  obtain  the  terms  in 
general  equation  (24).  Indeed,  when  a problem  containing  a 
particularized  example  is  selected,  the  expressions  of  these  are 
included  in  the  terms  of  (24). 

Current,  magnetic  flux,  electromotive  force  and  inductance 
definitions  have  been  obtained  using  some  terms  of  (24), 
simplifying  their  notation  and  speeding-up  their  computation 
during  the  solution  stage  rather  than  during  the  postprocessing 
stage. 

Finally,  the  approach,  implemented  in  a software  package 
developed  by  the  authors,  has  been  applied  to  several  cases. 
The  results  have  shown  a good  agreement  with  those 
measured  in  the  lab. 
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